A branch-point approximant for the equation of state of hard spheres 
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Using the first seven known virial coefficients and forcing it to possess two branch-point sin- 
gularities, a new equation of state for the hard-sphere fluid is proposed. This equation of state 
predicts accurate values of the higher virial coefficients, a radius of convergence smaller than the 
close-packing value, and it is as accurate as the rescaled virial expansion and better than the Pade 
[3/3] equations of state. Consequences regarding the convergence properties of the virial series and 
the use of similar equations of state for hard-core fluids in d dimensions are also pointed out. 



I. INTRODUCTION 

The virial expansion of the equation of state (EOS) is 
an expansion in powers of (usually) the number density 
p that was originally introduced phenomenologically by 
Kammcrlingh Onncs 1 in 1909 in order to provide a math- 
ematical representation of experimental data. Later, in 
what may be considered as one of the great achievements 
in statistical physics in the twentieth century, Mayer- 
was able to derive such an expansion for the pressure p 
of a classical fluid in terms of its density. The corre- 
sponding virial coefficients (usually denoted by Bj) turn 
out to be related to integrals over the interaction among 
groups of fluid particles and are in general functions of 
the absolute temperature T. In the case of hard-sphere 
(HS) fluids, which are the subject of this paper, the virial 
coefficients are however independent of T. In particu- 
lar, the value of the second virial coefficient for HSs in 
d dimensions is Bi = 2 d ~ 1 Vdcr d , where a is the diame- 
ter of the spheres and v d = (n/4) d / 2 /T(l + d/2) is the 
volume of a d-dimensional sphere of unit diameter, a re- 
sult first derived for three-dimensional HSs (d — 3) by 
van der Waals4 Analytical expressions for B3 and B4 
are also available in the Uteraturo 4 i 5 i 6 i 7 i 8 i 9 i 10 i u i 12 i 13 iM but 
higher virial coefficients must be computed numerically 
and, since this represents a non trivial task, up to now 
only values up to the tenth virial coefficient have been 
reported , 15 - 16 - 17 - 18 - 19 - 20 - 21 - 22 - 23 ' 24 - 25 - 26 - 27 ' 28 - 29 - 30 - 31 - 32 - 33 - 34 

The virial expansion for ci-dimensional HS systems is 
often expressed in terms of the packing fraction 7/ defined 
as 77 = Vd pa d . Hence, for these systems the compressibil- 
ity factor Z = p/ pksT (with ks the Boltzmann constant) 
is given by 



z(v) = 1 



bjrf 



3-1 



(1) 



where the (reduced) virial coefficients bj = Bj/(vd<J d y~ l 
are now pure numbers. 

The availability of only a few virial coefficients repre- 
sents a restriction on the usefulness of the virial expan- 
sion and many issues about it are still unresolved. For in- 
stance, its radius of convergence is not known eventhough 
lower bounds are availablej 35 ' 36 Secondly, although all the 
available virial coefficients in d = 2, d = 3, and d = 4 are 



positive, even the character of the series (either alternat- 
ing or not) is still unknown. In fact results from higher 
dimensions suggest that the positive character might not 
be true for the higher virial coefficients of hard disks and 
spheresi 33 i 37 ' 38 Finally, people have usually recurred to 
approximate EOSs obtained through the knowledge of 
the limited number of virial coefficients via various se- 
ries acceleration methods such as Pade or Levin approxi- 
mants. However, the expectation that these EOSs would 
ultimately lead to the complete phase behavior of the sys- 
tem has not been fulfilled. Hence, the question of whether 
the virial series contains relevant information related to 
the phase behavior of the HS system also remains as an 
open one. 

Recently it has been clearly established that the EOSs 
for hard hyperspheres (d > 4) predicted by the Percus- 
Yevick (PY) integral equation possess a branch-point sin- 
gularity on the negative real axis that is responsible for 
the radius of convergence and the alternating character 
of the virial seriesi 37 i 38 It is very likely that these features 
are not artifacts of the PY approximation but would be 
shared by the exact EOSs. However, in the case of hard 
spheres (d — 3), the radius of convergence of the PY 
EOS is artificially 77 = 1 and, as stated above, there is 
no definite indication about the nature of the singular- 
ity responsible for the true radius of convergence or its 
value M- 

The main aim of this paper is to shed some more 
light on the character of the virial series of the three- 
dimensional HS fluid. The idea is to propose a new 
(heuristic) EOS for HS systems in d dimensions that, 
for reasons that will become clear later, we will refer to 
as a 'branch-point approximant.' Such a proposal is not 
geared specifically towards obtaining an accurate EOS 
but rather relies on the notion that the radius of conver- 
gence of the virial series might be dictated by a branch- 
point singularity. In any case, the plausibility of this no- 
tion will be assessed by comparing the predictions of high 
virial coefficients coming out of the proposal both with 
the exact values of these coefficients for each d and with 
the performance of other proposals for the EOS (rescaled 
virial expansions and Pade approximants). 

The paper is organized as follows. In the next section 
we introduce the new EOS including a branch-point sin- 
gularity and examine the case of three-dimensional HSs. 
Section Hill refers to the use of the same type of EOS for 
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TABLE I: Expressions (for general k) and numerical values 
(for k = d = 3) of ai , <22, A, and ci— C3. 



Coefficient 



Expression 



Value 



A 

Cl 
C2 
C3 



7a? - 6Sk,7/Sk,s 



l(a 2 -ai) 2 /S KB 
3ai + 6 2 A 
§(a 2 + a?) - S fc , 3 4 
|ai(3a2 — a?) + Sk,iA 



0.27f232 
1.94804 
1.51486 
6.87314 
0.00268343 
1.33515 



different dimensionalities. We close the paper in Sect. II VI 
with some discussion and concluding remarks. 



II. THE CASE OF HARD SPHERES (d = 3) 

We begin by proposing a 'branch-point approximant' 
for the EOS of d-dimensional HS systems, namely 



1 + C177 + c 2 i) 2 + C377 3 - (l + 2ai?7 + a 2 rj 2 ) 
A(l - rf) k 



3/2 



z(v) = 1+ 

.n.i - // i 

(2) 

where A, a%, a 2 , c\, c 2 , C3, and k are parameters to be de- 
termined. This functional form (with k = 3) is inspired 
by the EOS for hard hyperspheres in d = 5 predicted 
by the PY theory through the virial route . 39 i 40 i 41 i 42 As 
stated above, we will assume the approximant form given 
in Eq. @ for three-dimensional HSs as a toy model to 
highlight the possibility that the radius of convergence 
of the virial series in this system might be dictated by a 
branch-point singularity. According to the philosophy of 
an approximant, the six coefficients A, a\, a 2 , ci, c 2 , and 
C3 are obtained from the knowledge of the virial coeffi- 
cients 62-67. The resulting expressions are given in Table 
HI where we have called 



E 

i=2 



k 

n-j 



(-iyb 3 



(3) 



Although the choice for k is in principle arbitrary, a 
natural one seems to take k — d. Hence, in this Sec- 
tion we assume k = 3. A special situation takes place if 
a 2 = a\ . In that case, the denominator (6*3,5) in the ex- 
pression for A must vanish in order to have a finite value, 
i.e., 65 = b 2 — 3&3 + 364. Since this denominator also ap- 
pears in the expressions for a\ and a 2 , the respective nu- 
merators (5*3,6 and 6*3,7) must also vanish, i.e., one must 
have be = 3b 2 — 863 + 664 and bj = 6b 2 — 1563 + IO64. Un- 
der those conditions, one has c\ = 3ai + b 2 A, c 2 = 3a? — 
5*3,3, A C3 = a\ + S3 ±A, so that Eq. © becomes Z(j]) = 
l + V [6a + (63 - 36 2 )77 + (64 - 36 3 + 36 2 )t7 2 ] /(l-vf, re- 
gardless of the values of a± and A. The aforementioned 
relationships are precisely satisfied by the virial and com- 
pressibility routes to the EOS in the PY approximation 



TABLE II: Exact and predicted values of bs-bw in the three- 
dimensional case 



Coefficient 


Exact 


Branch-point 
Eq. J5J, k = 3 


Rcscalcd expansion 
Eq. JS}, rn = 6, k = 3 


Pade [3/3] 
Eq. © 


bs 


68.538 


68.609 


68.812 


69.040 


69 


85.813 


85.532 


86.219 


87.147 


610 


105.78 


104.32 


105.56 


107.93 



for d — 3. Therefore, the functional form ([2]) is general 
enough as to include both PY EOSs, and thus also the 
Carnahan-Starling (CS) EOS^ given by 



Zos(v) 



1 + 77 + rf 



(4) 



as particular cases. Moreover, in the one-dimensional 
case one has bj = 1, so that again the relationships are 
satisfied and the resulting compressibility factor reduces 
to the exact EOS of the system, namely 2" (77) = 1/(1 — 7)). 

The numerical values of the coefficients ai, a 2 , A, 
and C1-C3 obtained from the known values of the first 
seven virial coefficients 30 ' 33 (namely, b 2 — 4, 63 = 10, 
64 ~ 18.364768, 65 ~ 28.2245, b 6 ~ 39.8151, b 7 ~ 
53.3444) are given in Table Q] The two branch points 
-(a 1 ±-s/aj-a 2 )/a 2 = -0.139234 ±0.702817?; lie on the 
complex plane. Their modulus is l/^/ai = 0.716 and this 
is then the radius of convergence of the virial series of the 
EOS @. While this radius is possibly an overestimate 
(in fact, it is larger than the freezing density), it is not un- 
physical since it is smaller than the close-packing value, 
in contrast to the radius 77 = 1 given by the PY, the 
CS, and the Carnahan-Starling-Kolafa (CSK)M EOSs, 
to name just a few. 

Table mi compares the know n 30 ' 33 and predicted values 
of 63-610 . Apart from the values predicted by Eq. $Z§, 
the table also includes the values obtained from the two 
following approximate EOSs that also make use of 62-67: 
the rescaled virial expansion^ 2 - 



1 + 

(I-77)* 



(5) 



k,n+l 



where m = 6, k — 3, and C n = (—1) 

with 5fe,„ given by Eq. ([3]), and the best^ Pade approx- 
imant [3/3] given by 



Z 



1 + ELiA.r?" 



(6) 



where D n and E n (n — 1,2,3) are combinations of 62-67 
whose explicit expressions may be easily obtained but will 
be omitted here. The deviations from the correct ones of 
the values for 6s-6io predicted by Eq. @ are 0.1%, 0.3%, 
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FIG. 1: (Color online) Plot of sgn(^) In \bj | for j < 150 in the 
three-dimensional case. 

and 1.4%, respectively. In contrast, the deviations of the 
values predicted by the rescaled virial expansion and the 
Pade approximant [3/3] are 0.4%, 0.5%, and 0.2% and 
0.7%, 1.5%, and 2%, respectively. Note that, in particu- 
lar, the rescaled virial expansion predicts a very accurate 
value for bio, even better than the prediction for bg. At 
a qualitative level, an interesting outcome of Eq. is 
first that it predicts a negative value of a certain coeffi- 
cient (specifically, 659) and secondly that henceforth the 
coefficients change sign every 1-2 terms. Figure [TJ shows 
sgn(foj) In \bj\ for j < 150. In contrast, the rescaled virial 
expansion predicts that all the bj are positive, while the 
Pade [3/3] predicts positive coefficients up to b$e and then 
alternating signs for groups of 55 consecutive coefficients. 

While the comparison between the exact values of the 
higher reduced virial coefficients and those that follow 
from the expansion of Eq. is quite satisfactory, one 
may reasonably wonder how the new EOS will perform 
when compared with other accurate proposals. Figure [2] 
shows that both the branch-point approximant and the 
rescaled virial expansion deviate less than 0.3% from the 
CS values for < 77 < 0.5 and are in very good agreement 
with simulation data 4 ^. The Pade [3/3] does a poorer 
job in this instance. Hence, the performance of the new 
proposal is also very accurate over the whole fluid phase 
range. 

III. OTHER DIMENSIONALITIES 

In this Section we perform a similar analysis of the use 
of Eq. © with k — d for all the dimensionalities (d = 
2, 4, 5, 6, 7, 8, and 9) where the first ten virial coefficients 



FIG. 2: (Color online) Plot of W0lZ(n)/Z C s (v) ~ 1] in the 
three-dimensional case. Solid line: Z{rf) given by the branch- 
point approximant ([2]) with k = 3; dotted line: Z(rj) given by 
the rescaled virial approximant <J75J with m = 6 and k — 3; 
dashed line: Z(n) given by the Pade [3/31 approximant JB). 
The circles are simulation data from Ref. liq . 



are known 1 33 1 3 

Table UTT] displays the exact and predicted values of 63- 
610 for d — 2 and d — 4-9 as given by the branch-point 
approximant, the rescaled virial expansion and the Pade 
[3/3] EOSs. 

For d = 2 the best approximant is the Pade [3/3]. The 
branch-point approximant also does a good job in the 
case of hard disks, but the rescaled virial approximation 
is slightly better (except for the value of 6g). The case 
d = 4 is somewhat peculiar because the predictions from 
all the approximants are rather poor. In any event, the 
Pade [3/3] gives the 'best' performance, followed by the 
rescaled virial expansion, and finally the branch-point 
approximant. This latter even 'anticipates' the likely 
aternating character of the series and predicts a nega- 
tive value of &10. The situation changes for 5 < d < 9 
where the performance of the rescaled virial expansion is 
extremely poor and in fact it never predicts negative coef- 
ficients, even when the exact be < is introduced (d = 6, 
d = 7) or the exact 64 < and be < are introduced 
(d = 8, d = 9). On the other hand, in these dimen- 
sionalities the Pade [3/3] predicts the right signs, while 
the branch-point approximant predicts, in addition, very 
good values. 

With regards to the EOS of hard disks, in Fig. [3] we 
compare the performance of the different approximants 
with respect to the simulation dat a 47 ' 48 in the range 
< i] < 0.7. Note that in this range the new proposal is 
able to capture the deviations (of up to 2%) of the simu- 
lation results that one gets from the use of the reasonably 
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TABLE III: 


Exact and predicted values ol os-oio 


for d = 2 and d = 4-9 




Coefficient 


Exact 


Branch-point 


Rescaled expansion 


Pade [3/3] 






Eq. ©, k = d 


Eq. m — 6, k — d 


Eq. © 






d = 2 






6s 


8.3191 


8.3397 


8.3408 


8.3241 


bg 


9.2721 


9.3711 


9.3297 


9.3001 


''hi 


10.216 


10.469 


10.319 


10.298 






d = 4 






bs 


605.66 


284.49 


486.07 


543.55 


bg 


739.88 


3339.4 


562.33 


605.51 


610 


1516.7 


-5388.0 


579.18 


704.91 






d = 5 






6s 


-3.0064 x 10 4 


-3.0177 x 10 4 


3.1662 x 10 4 


-2.6584 x 10 4 


bg 


3.2083 x 10 5 


3.1961 x 10 5 


9.4841 x 10 4 


2.3254 x 10 5 


610 


-3.3810 x 10 6 


-2.9014 x 10 6 


2.2311 x 10 5 


-1.6907 x 10 6 






d = 6 






6s 


-3.0752 x 10 7 


-3.0448 x 10 7 


9.4362 x 10 6 


-2.6002 x 10 7 


bg 


7.3370 x 10 s 


7.1562 x 10 s 


3.42934 x 10 7 


5.0061 x 10 8 


but 


-1.8472 x 10 10 


-1.7587 x 10 10 


9.3201 x 10 7 


-9.5796 x 10 9 






d = 7 






bs 


-8.7684 x 10 9 


-8.6759 x 10 9 


1.3044 x 10 9 


-7.6577 x 10 9 


bg 


4.7482 x 10 11 


4.6063 x 10 11 


5.3315 x 10 9 


3.3839 x 10 11 


610 


-2.7274 x 10 13 


-2.5651 x 10 13 


1.6168 x 10 10 


-1.4941 x 10 13 






d = 8 






6s 


-1.6114 x 10 12 


-1.5950 x 10 12 


1.2461 x 10 11 


-1.4369 x 10 12 


69 


1.8713 x 10 14 


1.8107 x 10 14 


5.6751 x 10 11 


1.3783 x 10 14 


&10 


—2.3160 x 10 16 


—2.1589 x 10 16 


1.9031 x 10 12 


— 1.3220 x 10 16 






d = 9 






6s 


-2.3219 x 10 14 


-2.2913 x 10 14 


9.6059 x 10 12 


-2.0919 x 10 14 


bg 


5.5879 x 10 1B 


5.3746 x 10 16 


4.8352 x 10 13 


4.1962 x 10 16 


bio 


-1.4436 x 10 19 


-1.3252 x 10 19 


1.7789 x 10 14 


-8.4195 x 10 18 



accurate EOS due to Henderson,— namely 

1 + ?y 2 /8 



(1 - nf 



(7) 



Concerning the nature of the singularities in these di- 
mensions, in Table Hvl we present the values of the singu- 
larity closest to the origin and the radius of convergence 
of the corresponding virial series, as predicted by the 
branch-point approximant ^ with k — d. For compar- 
ison, the radius predicted by the PY integral equation 
is also included in this table. One finds that the new 
proposal predicts complex branch points for 2 < d < 4. 
These are precisely the cases where all the known virial 
coefficients are positive. On the other hand, for 5 < d < 9 
the branch point closest to the origin is a negative real 
value. This agrees with the PY results, which gives some 
support to the alternating series scenario. Also note that 
the branch-point approximant and the PY radii of con- 
vergence tend to agree as d increases. 



TABLE IV: Singularity closest to the origin and radius of 
convergence of the virial series, as predicted by the branch- 
point approximant (J5J) with k = d. The radius predicted by 
the PY integral equation is also included. 



d 



Singularity 



Radius Radius (PY) 



2 0.3234 ± 0.4533i 0.557 1 

3 -0.1392 ±0.7028i 0.716 1 

4 -0.044223 ± 0.07526i 0.0873 0.15 

5 -0.07838 0.0784 0.057 

6 -0.02960 0.0296 0.024 

7 -0.01302 0.0130 0.011 

8 -0.006062 0.00606 0.0051 

9 -0.002925 0.00292 0.0024 



5 



- -0.5 



o 
o 




-1.0 



-1.5 



0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 

n 

FIG. 3: (Color online) Plot of W0[Z(r))/Z H (r)) - 1] in the 
two-dimensional case. Solid line: Z{rj) given by the branch- 
point approximant ([2]) with k = 2; dotted line: Z(rj) given by 
the rescaled virial approximant (0 with m — 6 and k — 2; 
dashed line: Z(t?) given by the Pade [3/3] approximant ([6]). 
The circles and squares are simulation data from Refs. [4?] and 
0, respectively. 

IV. CONCLUDING REMARKS 

In this paper we have introduced a new proposal for the 
EOS of a three-dimensional HS fluid which is built from 
the knowledge of the first seven virial coefficients and 
possesses two branch-point singularities in the complex 
plane. Although the choice we have made may appear 
to a certain extent arbitrary, it is perhaps the simplest 
one embodying the PY and CS EOSs for HSs in three 
dimensions as well as the exact Z in d = 1 and the PY 
virial EOS in d = 5. The same functional form was also 
assumed for the EOS of HS fluids in other dimensions. 
For d = 3 the new EOS predicts accurate values of the 
higher virial coefficients, a radius of convergence smaller 
than the close-packing value, and, irrespective of the fact 
that its construction did not aim at accuracy, it is very 
accurate when compared to simulation results and with 
other approximants involving the same number of known 
virial coefficients. This last feature was shown to be also 
shared by the two-dimensional case. The proposal is also 
robust with respect to small (~ 1%) deviations in the 
value of the seventh virial coefficient, certainly more ro- 
bust than either the rescaled virial expansion or the Pade 
[3/3]. 

Except for d = 4 (where, as already pointed out 
by Clisby and McCoy-2 3 . in a somewhat related con- 
text, perhaps one would require better accuracy of the 



known virial coefficients), in all other dimensionalities the 
branch-point approximant gives the best overall perfor- 
mance with respect to the prediction of the known virial 
coefficients. In particular, the rescaled virial expansion 
is unable to predict even the signs of known virial coef- 
ficients for d > 5 and the Pade [3/3], although correctly 
capturing these signs, leads to higher deviations. This of 
course constitutes no proof that the true EOS of HS sys- 
tems should include a branch-point singularity, but the 
evidence provided here is at least consistent with it. 

On a related vein, the new EOS for HSs in d = 3 also 
leads to an alternating virial series, with b§g being the 
first negative reduced virial coefficient. Given the dif- 
ficulty of computing exact high order virial coefficients, 
it is unlikely that the alternating series scenario for HSs 
in three dimensions may be confirmed in the near fu- 
ture. However, in view of the present results and those 
obtained in higher dimensions ; 33 ' 3 ' 38 it certainly gets re- 
inforced. 

One can reasonably wonder whether the present ap- 
proach to construct the EOS of HS systems using a num- 
ber of known virial coefficients may be cast in a sys- 
tematic way. While the answer is certainly not unique, 
the following constitutes a possible generalization. Wc 
rewrite the compressibility factor as 



1 + En=T c nV n - (1 + 2aiy + a-2V 
A(l - 7]) k 



2 ^Ar+l/2 



(8) 

where taking N — 1 corresponds to Eq. . In the three- 
dimensional case (A; = d = 3), the approximant with 
N = (which amounts to including only the first five 
virial coefficients) predicts be~bio with deviations equal 
to 0.54%, 0.69%, 1.2%, 3%, and 5%, respectively. On the 
other hand, the approximant with N — 2 predicts 610 
with a 1.2% deviation. Therefore, although there is cer- 
tainly an improvement in the prediction of 610 on going 
from N = to N = 2, a reasonable compromise between 
simplicity, generality, and accuracy seems to suggest that 
the choice N = 1 is the most adequate. 

Finally, it should be pointed out that if instead of 
choosing k = d as we have done in this paper, a different 
k is picked (say k — 3 for all d) we find slight variations in 
the numerical predictions but the overall picture remains 
unaltered. 
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